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PACS numbers: 71.10.-w, 02.70.Ss, 71.27.+a 

In a recent Letter Feldbacher, Held, and Assaad 
have proposed a new scheme for the Quantum Monte 
Carlo simulations of the impurity Anderson model. This 
method is supposed to calculate the average of an oper- 
ator O over the ground state of the Hamiltonian H via 
the expression 
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where Ht is an appropriate trial Hamiltonian (see RefQ, 
Eq.(3)). The main requirement to the choice of Ht is 
that its ground state can not be orthogonal to the ground 
state of the Hamiltonian under consideration H . In real 
calculations, the authors of Ref|jJ have taken H as the 
impurity Anderson Hamiltonian and as its nonin- 
teracting part, i.e. the Anderson Hamiltonian without 
the Hubbard interaction term —U/2 {n^ — where n„ 
is the operator of number of localized electrons with 
spin projection a =T)i on the impurity site. Unfortu- 
nately, this choice is not consistent with the basic restric- 
tion mentioned above due to the Anderson orthogonality 
catastrophe 2, 3]. Anderson has proven that in ther- 
modynamic limit iV — > oo (where N is the number of 
electrons) a local perturbation leads to complete recon- 
struction of the ground state of fermionic system in such 
a way that the overlap of the "old" and "new" ground- 
state way functions, |^'t) and j'l'), is proportional to 
N~°' where for the case of spherically symmetric pertur- 
bation potential 
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Si being the scattering phases at the Fermi level with 
orbital quantum number I. Later this expression was 
eneralized to the case of arbitrary local perturbation 


where S is the scattering matrix of the local perturba- 
tion at the Fermi surface. Original derivation of the An- 
derson orthogonality catastrophe was done for the case 
of local one-body potential; however, later it has been 
proven also for the Hubbard interaction term 0]. Tak- 
ing into account the Friedel sum rule for the Anderson 


model 1^ one can demonstrate that a is proportional 
to (fid ~ n^)^ where n^, are the impurity occupation 
numbers for the Hamiltonians H and Ht, correspond- 
ingly. For a special case of symmetric half-filled Ander- 
son model one has a = P|; however, to satisfy this in a 
generic case one has to know already an exact answer for 
which makes the whole scheme unpractical, at least, 
in the context of the dynamical mean-field theory. For 
the case of non-degenerate Anderson model the value 
can be calculated exactly by the Bethe Ansatz method 


Thus, for the infinite system the ground state does 
not contribute to both numerator and denominator of 
the ratio and, instead of being the average over the 
ground state the latter is rather some weighted average 
over excited states. 

The "orthogonality catastrophe" is a generic and fun- 
damental property of many-body fermionic systems (see, 
e.g., discussion in ReflT^ and should be carefully taken 
into account at the development of different projection 
techniques. 
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